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ABSTRACT 

Dirac sea corrections for bulk properties of finite nuclei are com- 
puted within a self-consistent scheme in the a-uo model. The valence 
part is treated in the Hartree approximation whereas the sea contribu- 
tion is evaluated semiclassically up to fourth order in K. Numerically, 
we find a quick convergence of the semiclassical expansion; the fourth 
order contributing much less than one percent to the binding energy 
per nucleon. 
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1 Introduction 



The quantum field theoretical approach to Nuclear Physics seems to be the 
only practical device known up to now to include relativistic corrections in the 
description of atomic nuclei in a systematic way [[I], |2|, [| [|, Q . In most works 
the description is limited to the valence part within a mean field Hartree 
approximation. Much less work has been devoted to the study of the Dirac 
sea of nucleons which should be included in principle to preserve the unitarity 
of the theory. The Dirac sea has been mostly considered in order to describe 
nuclear matter || or finite nuclei within the local density approximation 
up to two loops ||. On the other hand, even at the lowest order of the 
loop expansion, the physical relevance of the single particle negative energy 
nucleon states in finite nuclei is not fully understood. The explicit inclusion 
of the sea presents both numerical and conceptual problems. A numerical 
computation of the Dirac sea energy would require the diagonalization of 
all the continuum states of the single particle Hamiltonian and subsequent 
renormalization. A more serious problem is represented by the appearance 
of short distance tachyonic ghosts in the meson propagators at the one loop 
level which makes the mean field vacuum state |7], || and nuclear matter || 
unstable against formation of translationally non-invariant configurations of 
size of the order of 0.2 fm. In the relativistic approach to Nuclear Physics, a 
prescription has recently been suggested to deal with the instability problem 
based on the Kallen-Lehmann representation of two-point Green's functions 
which, so far, has only been applied to the study of nuclear matter and the 
corresponding local density approximation [|K], [|. In addition, the method 
presents ambiguities |! l]] and its extension to arbitrary Green's functions is 
not known. 

A way out of the above mentioned difficulties, which applies specifically 
to the mean field approximation, is to treat the Dirac sea assuming that the 
mesonic mean fields are smooth functions in coordinate space, thus allow- 
ing for a derivative or semiclassical expansion of the Dirac sea contribution 
to the energy, but keeping the shell structure in the valence part. Indeed, 
the problem of the high momentum instability of the sea is circumvented 
since only low Euclidean momenta are probed by the semiclassical expan- 
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sion. Moreover, the numerical problems related to the continuum states or, 
equivalently, the non-local nature of the fermionic determinant, are drasti- 
cally reduced to the study of an analytically given local functional of the 
mesonic fields and their lower derivatives. In fact, such an approach has 
been applied to compute the effect of the Dirac sea up to second order in 
h [0, [12| . There, it was found that such effects cannot be considered negli- 
gible. Obviously, there arises the question whether the higher order terms 
in the expansion can be safely neglected. There is, however, an additional 
subtlety, namely, even if the semiclassical expansion in powers of % turns out 
to be numerically convergent, there still could appear finite corrections of the 
type exp(— M^Rc/h) which might be numerically significant as compared to 
higher order h corrections. Such terms are beyond a semiclassical expansion 
and represent a sort of shell effects which, for the valence part, are known to 
be comparable to the second order h corrections Finally, an additional 



advantage of applying the semiclassical expansion to the Dirac sea lies in the 
absence of turning point divergences. In contrast, in the valence part these 
divergences take place and become a serious problem beyond second order, 
particularly in the formulation of a self-consistent approach. 

In the present letter, we undertake the calculation of the fourth order term 
of the Dirac sea contribution to the binding energy per nucleon as well as 
other relevant observables. This is done in a combined self-consistent Hartree 
treatment of the positive energy single particle states and a semiclassical 
expansion of the negative energy states. For simplicity, we do so for the 
o-uj model, although we expect our considerations to hold in more realistic 
versions including, e.g., p-meson as well as Coulomb interaction. 

2 Field theoretical model 

The a-u model is characterized by the following Lagrangian density [[![] 

C = ^[ 1 ^-g v V^)-(M-g s cft)}^+^(d^ ( ft-m 2 s( ft 2 ) + 

-\F,uF^ + l -mlV,V^ (1) 

where \1/ is the isospinor nucleon field, (ft the scalar field, the u;-meson field 
and F^ u = d^V u — d u V^. In the former expression the necessary counterterms 
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required by renormalization are implicitly included. In the time-independent 
mean field Hartree approximation the meson fields are replaced by classical 
c-numbers whereas the fermionic state is given by a Slater determinant con- 
sisting of valence plus sea single particle states. The total mean field energy 
is given by E[<f>, V ] = E val + £ sea + E B , where 



n 

"(V</>) 2 + m s( j) 2 - (VV ) 2 - m v V 2 } ■ (2) 



E B = I I d 3 x 



2 

and the single particle valence and sea orbitals depend functionally on the 
mesonic fields <fi and Vq through the Dirac equation 

[-ia-V + g v V (r)+(3(M-g s (f>(r))}i; n (r) = E n il> n (r)- (3) 

Notice that in the nuclear ground state of spherical nuclei, the only case to be 
considered in this letter, the spatial components of the w-meson field vanish 

1- 

To perform in practice the derivative expansion we borrow from refs. |TJ 



13| , where the single particle level density p(E) is obtained to fourth order in 
h for an arbitrary space dimension D. The sea energy can be deduced from 
eqs. (3) and (18) of ref. [[14] by means of the relation 



£ sea = f dEEp(E). (4) 

The sea energy is ultraviolet divergent and requires a suitable renormaliza- 
tion. The standard dimensional regularization can be applied with D = 
3 — 2e. Since we have in mind a semiclassical expansion of the sea energy we 
specify the renormalization conditions at zero momentum as done in ref. [||]. 
This results in the following expression for the renormalized sea energy to 
fourth order in the gradients of the mesonic fields 



ITisca 
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£ "* = i^/ dSr {5 io 4< w > 2 - io 4< v *> 2 ) 

ET a = 7 - / d 3 r { - 11 $- 4 (V$) 4 - 22 $- 4 (W) 2 (V$) 2 
57607T 2 J [ 

+44 $- 4 ((Vi$)(V^)) 2 - 44 $- 3 ((V i $)(V i V))(V 2 ^) 
-8$" 4 (VV) 4 + 22$- 3 (V 2 $)(V$) 2 + 14$- 3 (VV) 2 (V 2 $) 

-18$~ 2 (V 2 $) 2 + 24$~ 2 (VV) 2 j . (5) 

Here, V A = g„Vo,$ = M — g s (p and 7 is the spin and isospin degeneracy 
of the nucleon, i.e., four if there are two nucleon species. The first two 
terms have been known for some time, see e.g. [EL The fourth order term 
is new and turns out to be ultraviolet convergent. This was expected in 
advance since dimensional counting shows that there is a correspondence 
between the number of gradients and the superficial degree of divergence of 
the corresponding momentum integral. 



3 Numerical results 

The mean field equations of motion are obtained, as usual, by minimizing 
the total mean field energy functional, E[<fi, Vq\ (see definition above eq. (0)), 
with respect to arbitrary variations of 0(r) and Vb(r). This yields eq. @ 
together with 

(V 2 -m 2 )0(r) = -9s(pT\r) + pT(r)) 
(V 2 -ml)V (r) = -^(p-i( r ) + P r(r)) (6) 

where we have explicitly separated the valence and sea contributions to the 
scalar and baryonic densities, p s = and p v = (^^ r ), respectively. For 

both densities, the valence part is just the corresponding sum over valence 
orbitals. On the other hand, the sea densities can be obtained as 

1 5E sea v . . 1 5E sea 

9s <w) g v SV {r) 

In turn, the expansion for .E sea given in eq. @ provides a similar semiclassical 
expansion for the scalar and baryonic densities. We will consider closed-shell 
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nuclei for which the equations admit spherically symmetric solutions. These 
equations must be solved self-consistently. This poses no particular problems 
when the sea is included to zeroth or second orders since they give rise to 
second order differential equations. The fourth order case is qualitatively 
different since it requires to solve fourth order differential equations and 
consequently further boundary conditions need to be specified. Instead of 
doing so, the fourth order sea contributions are treated perturbatively, that 
is, we consider the full self-consistent solution to second order and use the 
resulting meson fields to compute the fourth order sea corrections to the 
binding energy and the nuclear radius. Obviously, due to the stationarity of 
the energy functional to second order, this procedure reproduces the correct 
fourth order correction to the energy up to higher order terms. As it will be 
shown below, the fourth order corrections turn out to be small, thus justifying 
a posteriori this procedure. 

In Table [I] we present results for the binding energy per nucleon and mean 
squared charge radius for several closed-shell nuclei calculated in various ap- 
proximations. In all cases, the parameters of the Lagrangian, namely, g s , g v , 
m s and m v , are adjusted to reproduce the u- meson mass, the saturation den- 
sity and binding energy per nucleon for nuclear matter and the mean squared 
charge radius of 40 Ca [@, [5|. The corresponding numerical values are given 
in Table |2|. In both tables the entry "no sea" stands for the customary va- 
lence Hartree approximation, where the effects due to the Dirac sea are fully 
neglected. Oth- and 2nd-order represent the results of successive inclusion of 
the Dirac sea keeping the corresponding orders in the expression for the sea 
energy and sea densities in eqs. fl5|) and (0). The experimental values are 
given only for illustration, since the a-u model must be supplemented with 
additional degrees of freedom, such as the p-meson and Coulomb interaction, 
to be realistic || ||. From Table [I], one can see that the Dirac sea corrections 
are not globally small and hence cannot be neglected. The effect becomes 
more dramatic on the fixing of the parameters, Table Finally, in Table |3] 
we present our results for the full calculation where the Dirac sea is described 
semiclassically up to fourth order in h. The various contributions to minus 
the total binding energy per nucleon are displayed, namely, valence kinetic 
energy, valence potential energy, minus sea potential energy and total sea 
energy. Their total sum yields minus the binding energy per nucleon. 

In Table p], the results for the fourth order have been intentionally omitted 
since they do not further modify the presented second order results. In fact, in 
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all cases considered, i.e., closed-shell nuclei, the fourth order contributions are 
very small as they do not significantly influence neither the bulk properties 
nor the parameters of the model. The smallness of the correction does not 
follow from big cancelations of any kind; each term in eq. ([|) turns out to 
be uniformly small in the integration region. This is an important result 
for it indicates that the semiclassical expansion converges much faster than 
one might naively expect. Indeed, the fourth order sea correction to the 
binding energy per nucleon can be estimated, from eq. (|5]), to be of the 
order of 1/(RA), R being the nuclear radius and A the mass number. In 
the case of calcium it corresponds to a correction of about 1 MeV. The 
fact that this number turns out to be much smaller is a direct consequence 
of the overall numerical dimensionless coefficient, whose value requires an 
explicit calculation, as done here. The higher orders are not expected to be 
relevant. For instance, the sixth order can be estimated to be suppressed by 
a further factor 1/(U R) 2 , U being the depth of either the scalar or the vector 
potentials, i.e., about one percent of the fourth order correction for calcium. 

As already mentioned in the introduction, a word of caution is in or- 
der. The fast numerical convergence of the semiclassical expansion for the 
sea contribution to the mean field energy does not necessarily imply that it 
converges to the exact Hartree result. The expansion might be only asymp- 
totic or it might be convergent but not to the exact result. Quite generally, 
non analytical terms in h, which are beyond a semiclassical approach, are 
expected to occur in global properties. Actually, such contributions can be 



shown to appear in simple non relativistic quantum mechanical models [15 



The existence of these terms must be kept in mind if the semiclassical results 
are to be compared with an exact calculation, since they might introduce 
systematic deviations. Such deviations seem to have been observed in the 



quark chiral soliton model [16 



4 Conclusions 

We summarize our points. The Dirac sea corrections to the mass and size of 
finite nuclei have been computed in a semiclassical self-consistent treatment 
of the a-oj model, but taking into account the shell structure of the valence 
part by means of a Hartree approximation. In this combined self-consistent 
approach, we have investigated the numerical convergence of the semiclas- 



7 



sical expansion up to fourth order in h for the considered bulk properties. 
Although the Dirac whole cannot be considered negligible, the over- 

all effect becomes smaller after the necessary readjustment of the parameters 
in the model. Furthermore, a fair description of the sea can be achieved 
by including just the zeroth and second order semiclassical terms. We em- 
phasize the non triviality of this result since a simple dimensional estimate 
allows much larger corrections than the ones actually found at fourth order. 
Finally, the present calculations can be regarded as a previous step to an 
exact mean field description of the polarization of the Dirac sea. 
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Table Captions 



1. Binding energy per nucleon (in MeV) and mean squared charge radius 
(in fm) for closed-shell nuclei in successive calculational schemes. The 
"no sea" column stands for the purely valence calculation. The "sea 
Oth" and "sea 2nd" columns refer to inclusion of the sea through a semi- 
classical expansion. In all cases the parameters have been readjusted 
as explained in the main text. Experimental values are only given as 
orientative reference. 

2. Parameters of the a-u model adjusted to fit the mean squared charge 
radius of 40 Ca and nuclear matter saturation density and binding en- 
ergy (see ref. || for the definition of the parameters). The nucleon and 
oj meson masses are fixed to 939 MeV and 783 MeV respectively. The 
same calculational schemes as in the Table [I] are considered. The scalar 
meson mass is given in MeV. 

3. In different columns, the valence kinetic energy, the valence potential 
energy, minus the sea potential energy and the total sea energy per 
nucleon are displayed for the full fourth order semiclassical approxima- 
tion to the Dirac sea. Their total sum yields minus the binding energy 
per nucleon. 
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No sea 


Sea 0th 


Sea 2nd 


Exp. 




B/A 


m.s.c.r. 


B/A 


m.s.c.r. 


B/A 


m.s.c.r. 


B/A 


m.s.c.r. 


20 Ca 


6.28 


3.48* 


6.00 


3.48* 


6.33 


3.48* 


8.55 


3.48 


48 p 


6.52 


3.47 


6.03 


3.51 


6.34 


3.52 


8.67 


3.47 


INi 


7.24 


3.72 


6.51 


3.79 


6.80 


3.79 


8.64 




90 ry 


8.36 


4.22 


7.99 


4.23 


8.22 


4.25 


8.71 


4.27 


^ 2 Sn 


8.81 


4.60 


8.43 


4.66 


8.62 


4.68 


8.36 




2 iPb 


9.84 


5.35 


9.55 


5.39 


9.70 


5.41 


7.87 


5.50 



Table 







m s 


9s 


cl 


9v 


Nosea 
SeaOth 
Sea2nd 


357.741 
227.840 
227.840 


449.93 
368.37 
406.6 


9.06283 
5.92153 
6.53607 


274.106 
147.524 
147.524 


13.8056 
10.1281 
10.1281 



Table g 





K val 


-|p val U 


l p sea u 


rpsea 
^0 


gsea 




40 p„ 


14.6 


-18.0 


-5.29 


1.68 


0.681 


-2.6-10" 


-4 


48 p 


15.6 


-18.8 


-5.61 


1.82 


0.670 


8.3-10" 


-5 


INi 


16.1 


-19.5 


-6.11 


2.01 


0.679 


6.7-10" 


-4 


l> 


16.1 


-20.3 


-7.02 


2.42 


0.586 


-1.4-10- 


-5 




16.7 


-21.0 


-7.52 


2.66 


0.526 


-l.o-io- 


-4 


2 iPb 


16.8 


-21.7 


-8.31 


3.01 


0.455 


-1.9-10 


-4 



Table | 
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